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Abstract 

The shape of jets produced in quasi-real photon-proton collisions at centre-of-mass 
energies in the range 134 — 277 GeV has been measured using the hadronic energy flow. 
The measurement was done with the ZEUS detector at HERA. Jets are identified using 
a cone algorithm in the tj — <j> plane with a cone radius of one unit. Measured jet shapes 
both in inclusive jet and dijet production with transverse energies E J T et > 14 GeV arc 
presented. The jet shape broadens as the jet pseudorapidity (rf et ) increases and narrows as 
increases. In dijet photoproduction, the jet shapes have been measured separately for 
samples dominated by resolved and by direct processes. Leading- logarithm parton-shower 
Monte Carlo calculations of resolved and direct processes describe well the measured jet 
shapes except for the inclusive production of jets with high if et and low E™ . The observed 
broadening of the jet shape as rf et increases is consistent with the predicted increase in 
the fraction of final state gluon jets. 
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1 Introduction 



Photoproduction at HERA is studied via ep scattering at low four-momentum transfers (Q 2 ~ 0, 
where Q 2 is the virtuality of the exchanged photon). In photon-proton reactions, two types 
of QCD processes contribute to jet production at leading order (LO) |], : either the photon 
interacts directly with a parton in the proton (the direct process) or the photon acts as a 
source of partons which scatter off those in the proton (the resolved process). The first year 
of HERA operation led to the observation of hard scattering in jp collisions with evidence for 
multijet structure as well as the presence of the resolved process 0, ||. Measurements of dijet 
events allowed the separation of the resolved and direct processes ||. The jet profiles have been 
measured and found to be described by leading-logarithm parton-shower calculations, except 
for resolved processes in the forward region || |7], |fl. 

In this paper, the internal structure of photoproduced jets is studied at the hadron level. 
The investigation of the internal structure of jets gives insight into the transition between a 
parton produced in a hard process and the experimentally observable spray of hadrons. In the 
present study, jets are searched for with an iterative cone algorithm H [TOj] with radius R = 1 
in the pseudorapidityQ (rj) - azimuth (</>) plane. The jet shape is defined as the average fraction 
of the jet's transverse energy (E^) that lies inside an inner cone of radius r concentric with 
the jet defining cone JTTJ] : 

^(r) = ^V E ^ r) , (1) 
1 ' N jets fe s E T (r = R) 1 1 

where Et{t) is the transverse energy within the inner cone of radius r and Nj ets is the total 
number of jets in the sample. By definition, ip(r — R) — 1. 

The jet shape is affected by fragmentation and gluon radiation. However, at sufficiently 
high Ej 1 the most important contribution is predicted to come from gluon emission off the 
primary parton and, therefore, is calculable in perturbative QCD. The lowest-non-trivial order 
contribution to the jet shape is given by next-to-leading order (NLO) QCD calculations for the 
reaction AB — > jet + X. Perturbative QCD predicts that gluon jets are broader than quark 
jets as a consequence of the fact that the gluon-gluon is larger than the quark-gluon coupling 
strength. 

Measurements of jet shapes have been made in pp collisions at y^i = 1.8 TeV using only 
charged particles |T^] as well as both neutral and charged particles |T3| , and a qualitative 



agreement with NLO QCD calculations [11, 14] was found. Similar measurements have been 



made in e + e~ interactions at LEP1 using both neutral and charged particles |TS| and found 
to be well described by leading-logarithm parton-shower Monte Carlo calculations. It was 
observed JTH] that the jets in e + e~ are significantly narrower than those in pp. This is due to 
the different mixtures of quark and gluon jets in these two environments [15] . Measurements 



of the jet width at LEP1 have shown that gluon jets are indeed broader than quark jets H . 

In this paper, measurements are presented of the jet shapes in both inclusive jet and dijet 
photoproduction at centre-of-mass energies in the range 134 — 277 GeV. The data sample used 
in this analysis was collected with the ZEUS detector in e + p interactions at the HERA collider 
and corresponds to an integrated luminosity of 2.65 pb _1 . Jets are selected with E%T > 14 GeV 
and — 1 < i]^ et < 2. The jet shape is measured using the ZEUS calorimeter and corrected to the 



x The ZEUS coordinate system is defined as right-handed with the Z-axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X-axis horizontal, pointing towards the centre of HERA. 
The pseudorapidity is defined as r\ = — ln(tan | ) , where the polar angle 9 is taken with respect to the proton 
beam direction. 
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hadron level. The measurements are presented as a function of the jet transverse energy and 
pseudorapidity. In dijet photoproduction, the jet shapes are measured for resolved and direct 
processes defined by the fraction of the photon's momentum participating in the production of 
the two jets of highest Ej? . 

Measurements of jet shapes are compared to QCD calculations based on different approaches: 

1) LO QCD calculation including initial and final state QCD radiation in the leading- 
logarithm parton-shower approximation as implemented in the program PYTHIA 5.7 



r7|. The final state parton system is hadronised using the LUND string model [IS 



LO QCD calculation as implemented in the program PYTHIA without initial and final 
parton radiation. The final state parton system is hadronised using the LUND string 
model. 



3) Fixed-order perturbative QCD calculation of the reaction e + p — ► 3 partons + X fl9| , |20 
Fragmentation effects are not included. 

These comparisons allow the study of the relative importance of parton radiation and frag- 
mentation in the formation of a jet as well as the differences between quark and gluon jets. It 
should be noted that the first two predictions refer to the hadron level while the third refers to 
the parton level. For the first two predictions, jets are searched for in the final state hadronic 
system using the same jet algorithm as in the data. For the third prediction, the experimental 
jet algorithm has been simulated by the introduction of an additional parameter, Rsep, as 



proposed in [I I 



2 Experimental conditions 

During 1994 HERA operated with 153 colliding bunches of protons of energy E p = 820 GeV 
and positrons of energy E e = 27.5 GeV, with 96 ns between bunch crossings. 



A description of the ZEUS detector can be found in |2T|, ^2fl. The components used in this 
analysis are briefly discussed. The uranium-scintillator calorimeter (CAL) E3] covers 99.7% of 
the total solid angle. It consists of the forward calorimeter (FCAL) covering the range 2.6° < 
9 < 36.7° in polar angle, the barrel calorimeter (BCAL) covering 36.7° < 9 < 129.1°, and 
the rear calorimeter (RCAL) covering 129.1° < 9 < 176.2°. Holes of 20 x 20 cm 2 in the centre 
of FCAL and RCAL are required to accommodate the HERA beampipe. For normal incidence, 
the depth of the CAL is 7 interaction lengths in FCAL, 5 in BCAL and 4 in RCAL. Each of 
the calorimeter parts is subdivided into towers which in turn are segmented longitudinally into 
one electromagnetic (EMC) and one (RCAL) or two (FCAL, BCAL) hadronic (HAC) sections. 
The sections are further subdivided into cells of inner-face sizes of approximately 5 x 20 cm 2 
(10 x 20 cm 2 in the RCAL) for the EMC and 20 x 20 cm 2 for the HAC sections. Each cell 
is viewed by two photomultipliers. At 9 = 90° the size of an EMC (HAC) cell in the f] — <p 
plane is approximately 0.04 x 11° (0.16 x 11°). Under test beam conditions the calorimeter has 



an energy resolution of a/E = 18%/yE(GeV) for electrons and ojE = 35%/jE(GeV) for 
hadrons. In order to minimise the effects of noise due to the uranium radioactivity, all EMC 
(HAC) cells with an energy deposit of less than 60 (110) MeV are discarded in the analysis. For 
isolated energy deposits, consisting of one cell surrounded by empty cells, this cut was increased 
to 100 (150) MeV. Particles impinging on the CAL lose energy in the inactive material in front 
of the CAL. In the region relevant for the present analysis, the inactive material constitutes 
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about one radiation length except in the region around the rear beampipe, 9 > 170°, and the 
solenoid support structure, 25° < 9 < 45° and 130° < 9 < 145°, where it is up to 2.5 radiation 
lengths. For the following measurements, the transverse energy and the shape of the jets have 
been corrected for these energy losses (see Section 5). 

The tracking system consists of a vertex detector (VXD) p4"| , a central tracking chamber 
(CTD) and a rear tracking detector (RTD) |22[ enclosed in a 1.43 T solenoidal magnetic 
field. The interaction vertex is measured with a typical resolution along (transverse to) the 
beam direction of 0.4 (0.1) cm. 

Proton-gas events occurring upstream of the nominal interaction point are out of time with 
respect to the e + p interactions and are rejected by timing measurements using a set of scintil- 
lation counters. 



2.1 Trigger conditions 



The ZEUS detector uses a three-level trigger system [^] . At the first level events were triggered 
on a coincidence of a regional or transverse energy sum in the CAL and at least one track from 
the interaction point measured in the CTD. At the second level a total transverse energy of 
at least 8 GeV, excluding the energy in the eight CAL towers immediately surrounding the 
forward beampipe, was required, and cuts on CAL energies and timing were used to suppress 
events caused by interactions between the proton beam and residual gas in the beampipe. 

The full event information was available at the third- level trigger (TLT). Tighter timing 
cuts as well as algorithms to remove beam-halo muons and cosmic muons were applied. For 
this analysis, the following additional conditions were required: a) the event has a vertex 
reconstructed by the tracking chambers with the Z value in the range \Z\ < 60 cm; b) E—pz > 
8 GeV, where E is the total energy as measured by the CAL, E = J2iEi, and pz is the Z- 
component of the vector p = Y^i E%^%\ i n both cases the sum runs over all CAL cells, Ei is the 
energy of the calorimeter cell i and fj is a unit vector along the line joining the reconstructed 
vertex and the geometric centre of the cell i; c) Pz/E < 0.95 to reject beam-gas interactions 
(this cut was not applied for events with E —p z > 12 GeV); and d) the total transverse energy 
as measured by the CAL, excluding the cells with polar angles below 10°, exceeds 20 GeV. 

For studies of the trigger efficiency an additional sample of events was selected by the TLT 
based upon jets found using a cone algorithm with radius R = 1 applied to the CAL cell energies 
and positions. The events were required to fulfill the same conditions a), b) and c) as above, 
and to have at least one jet of transverse energy, as measured by the CAL, E^l^ > 6.5 GeV 

and Vial < 2-5. 



3 Data selection 

Events from quasi-real photon proton collisions containing jets were selected offline using similar 
criteria as reported previously ||. The main steps are briefly discussed here. A search for jet 
structure using the CAL cells (see Section 5) is performed, and events with at least one jet of 
transverse energy, as measured by the CAL, E^lai > 10 GeV and —1 < rf^ < 2 are retained. 
The contamination from beam-gas interactions, cosmic showers and halo muons is negligible 
after demanding: a) at least two tracks pointing to the vertex; b) the vertex position along the 
beam axis to lie in the range —29 cm < Z < 36 cm; c) less than five tracks not associated with 
the vertex and compatible with an interaction upstream in the direction of the proton beam; 
d) the number of tracks not associated with the vertex be less than 10% of the total number 
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of tracks; and e) the total missing transverse momentum (j/t) be small compared to the total 



transverse energy (E^) by requiring j/t/JE^ 1 < 2 GeVa. Deep inelastic e + p scattering (DIS) 
neutral current events with an identified scattered positron candidate in the CAL according to 
the algorithm described in || are removed from the sample. 

The selected sample consists of events from e + p interactions with Q 2 < 4 GeV 2 and a median 
of Q 2 ~ 10~ 3 GeV 2 . The jp centre-of-mass energy (W) is calculated using the expression 
W = a/^s, where y is the inelasticity variable and s is the squared e + p centre-of-mass energy 
(300 2 GeV 2 ). The event sample is restricted to the kinematic range 0.2 < y < 0.85 using the 



following procedure. The method of Jacquet-Blondel fl26| , y,jB = (E — pz)/(2E e ) } is used to 
estimate y from the energies measured in the CAL cells. Due to the energy lost in the inactive 
material in front of the CAL and to particles lost in the rear beampipe, yjB systematically 
underestimates the true y by approximately 20%, an effect which is adequately reproduced 
in the Monte Carlo simulation of the detector. To compensate for this deficiency, the event 
selection required 0.16 < y JB < 0.7. The data sample consists of 15,368 events with a total of 
18,897 jets. The only significant remaining background is from unidentified DIS neutral current 
interactions with Q 2 > 4 GeV 2 , which is estimated by Monte Carlo techniques to be below 2%. 

4 Monte Carlo simulation 

The response of the detector to jets and the correction factors for the jet shapes were determined 
from samples of Monte Carlo events. 



The programs PYTHIA 5.7 |17[ and HERWIG 5.8 |27] were used to generate photoproduction 
events for resolved and direct processes. In PYTHIA the positron-photon vertex was modelled 
according to the Weizsacker- Williams approximation. In the case of HERWIG, the exact matrix 
elements were used for direct processes (e + g — > e + qq and e + q — > e + qg) and the equivalent 



photon approximation for resolved processes. Events were generated using GRV-HO [28] for 



the photon parton distributions and MRSA |29| for the proton parton distributions. In both 



generators, the partonic processes were simulated using LO matrix elements, with the inclusion 
of initial and final state parton showers. Fragmentation into hadrons was performed using 



the LUND string model |T8| as implemented in JETSET [^] in the case of PYTHIA, and the 
cluster model in the case of HERWIG. Samples of events were generated with different values of 
the cutoff on the transverse momentum of the two outgoing partons, starting at PTmin = 8 GeV. 

Additional samples of events were generated using the option of multiparton interactions 
(MI) in PYTHIA. This option, which applies only to resolved processes, adds interactions 
between the partons in the proton and the photon remnants calculated as LO QCD processes 
to the hardest scattering process of the event. The PYTHIA MI events were generated with a 
cutoff for the effective minimum transverse momentum for multiparton interactions of 1.0 GeV 
and with a cutoff on the transverse momentum of the two outgoing partons from the hardest 
scattering of pTmin = 8 GeV. 

All generated events were passed through the ZEUS detector and trigger simulation programs 



221. They were reconstructed and analysed by the same program chain as the data. 



5 Jet search and reconstruction of the jet shape 

An iterative cone algorithm in the rj — <fi plane |)], [K]] (PUCELL) is used to reconstruct jets, 
from the energy measured in the CAL cells for both data and simulated events, and also from 
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the final state hadrons for simulated events. 

The procedure is explained in detail for the jet reconstruction from the CAL cell energies {cal 
jets). In a first step, each CAL cell with a transverse energy in excess of 300 MeV is considered 
as a seed for the search. Their corresponding rj — <fi values are obtained from the unit vectors 
joining the vertex of the interaction and the geometric centres of the cells. The seeds are then 
combined into preclusters starting from that with highest transverse energy if their distance in 



the n — cj) plane, J (An) 2 + (A0) 2 , is smaller than 1 unit. The axis of the precluster is defined 
according to the Snowmass convention fll0"|] , where rjP reduster ^reciuster^ - g ^ e transverse-energy 
weighted mean pseudorapidity (azimuth) of all the seeds belonging to that precluster. 

In a second step, a cone of radius R = 1 is drawn around each precluster and all the CAL cells 
within that cone are combined to form a cluster. The axis of the cluster is defined according to 
the same prescription as for the preclusters but including all the CAL cells belonging to that 
cluster. A new cone of R = 1 is then drawn around the axis of the resulting cluster. All cells 
with geometric centres inside the cone are used to recalculate a new cluster axis. The procedure 
is iterated until the content of the cluster remains unchanged. 

In a third step, the energy sharing of overlapping clusters is considered. Two clusters are 
merged if the overlapping energy exceeds 75% of the total energy of the cluster with the lower 
energy; otherwise two different clusters are formed and the common cells are assigned to the 
nearest cluster. Finally, a cluster is called a jet if the transverse energy as measured by the 
CAL, -E^cah exceeds 10 GeV. The angular variables associated with the cal jets are denoted 
by vlai and 

The following procedure was used to reconstruct the jet shape from the CAL cells: for each 
jet the sum of the transverse energies of the CAL cells assigned to the jet with a distance 



(Arj) 2 + (A(f)) 2 to the jet axis smaller than r is determined, ET, C ai(r) } and divided by 
Et,co,i(j' = 1)- The jet shape as measured with the CAL, ip ca i(r), is thus defined as: 

where the sum runs over all the jets in the selected sample and Nj ets is the total number of jets 
in the sample. 

For the Monte Carlo events, the same jet algorithm is also applied to the final state particles. 
In this search, all particles with lifetimes longer than 10~ 13 s and with polar angles between 
5° and 175° are considered. The jets found are called hadron jets and the variables associated 
with them are denoted by E^ hadl rf^ ad , and 4^ d - Hadron jets with E^ had > 14 GeV and 
— 1 < n^ d < 2 are selected. The same jet shape definition as used above for the CAL cells is 
applied to the final state particles in the case of simulated events and the resulting jet shape is 
denoted by i/>f£[ (r). 

The comparison of the reconstructed jet variables between the hadron and the cal jets in 
simulated events [0] shows no significant systematic shift in the angular variables rf^ and </>^ 
with respect to n J h e ad and (pj^d- However, the transverse energy of the cal jet underestimates that 
of the hadron jet by an average amount of ~ 16% with an r.m.s. of 11%. The transverse energy 
corrections to cal jets averaged over the azimuthal angle were determined using the Monte Carlo 
events ||. These corrections are constructed as multiplicative factors, CiE^^n 3 ^), which, 
when applied to the Et of the cal jets give the 'corrected' transverse energies of the jets, 
£^f* = C(E^ t cal ,n 3 ca i l ) x E^l al ||. These corrections mainly take into account the energy losses 
due to the inactive material in front of the CAL. 
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5.1 Jet shape corrections 

The jet shapes as measured with the CAL are corrected to the hadron level using the Monte 
Carlo event samples. The corrected jet shapes are denoted by ip(r) and refer to jets at the 
hadron level with a cone radius of one unit in the t] — <p plane. The measurements are given 
for jets of corrected transverse energy E J T et > 14 GeV and -1 < nf et < 2, and in the kinematic 
region defined by Q 2 < 4 GeV 2 and 134 < W < 277 GeV. 

The reconstructed jet shapes are corrected for acceptance and smearing effects using the 
samples of Monte Carlo events of resolved and direct processes. The correction factors also 
take into account the efficiency of the trigger, the selection criteria, the purity and efficiency 
of the jet reconstruction, and the effects of energy losses due to the inactive material in front 
of the CAL. The corrected jet shapes are determined bin-by-bin as ip(r) = F^f (r) ■ ipcaiij'), 
where the correction factors F^f(r) are defined as F^ c (r) = ^f^(r) /^f (r). F^f(r) is 
determined separately for each interval in rf et and E 3 ^ 1 . 

For this approach to be valid, the uncorrected jet shapes in the data must be described 
by the Monte Carlo simulations at the detector level. As shown later, this condition is in 
general satisfied by both the PYTH1A and HERWIG simulations although some disagreement 
is observed in the forward region in the case of inclusive jet production. The latter can be 
reduced by adjusting the relative contribution of direct processes and the fraction of gluon jets 
in both resolved and direct processes. In the simulated events a cat jet is classified as a quark 
or gluon jet depending on the type (quark or gluon) of the closest parton from the two-to-two 
hard subprocess. 

The following procedure was adopted to obtain the best description of the uncorrected jet 
shapes: a) the relative contributions of resolved and direct processes were tuned by a fit to the 
measured RCAL energy distribution in the data; a distinct distribution of the energy deposit 
in the rear direction for resolved (direct) processes is expected due to the presence (absence) of 
the photon remnant; b) the fraction of gluon jets in direct and resolved processes was adjusted 
by a fit to the uncorrected jet shape. This procedure was applied separately for each rf et and 
E^ et interval. The correction factors F™f(r) are taken from this tuned version of PYTHIA. 
For the comparisons to the measurements presented in Sections 6 and 7 the untuned version of 
PYTHIA has been used. 

The correction factors do not show a strong dependence on rf et or Ej? and differ from unity 
by less than 10% (5%) for r > 0.4 (r > 0.6). For r = 0.3 (r = 0.2) they can increase to 
15% (30%). The dependence of the correction factors on the choice of fragmentation scheme, 
relative contributions of direct and resolved processes, fractions of gluon and quark jets in each 
process, and the inclusion of multiparton interactions in resolved processes was studied and 
found to be small. The resulting differences are taken into account as contributions to the total 
systematic uncertainty assigned to the measurements reported in the next sections. 

5.2 Systematic uncertainties of the measurements 

A detailed study of the sources contributing to the systematic uncertainties of the measurements 
was carried out. For each source, the largest change in the corrected value of ip(r) (indicated in 
parentheses) at a fixed value of r typically occurs at r = 0.2 and the magnitude of the induced 
change decreases rapidly as r increases: 

• The correction functions to the jet shapes were calculated using the event samples of the 
untuned version of PYTHIA (± 3%). 
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• Using the HERWIG generator to evaluate the energy corrections to the jets and the 
correction functions to the jet shapes (± 5%). 

• The correction functions were calculated using Monte Carlo samples of PYTHIA events 
consisting exclusively of either gluon or quark jets (± 9%). 

• The use of the PYTHIA generator including multiparton interactions in resolved processes 
(±4%). 

• Variation of the absolute energy scale of the CAL in the simulated events by ± 3% (± 2%). 

• The uncertainty in the simulation of the CAL response to low-energy particles (—2%). 

• Uncertainties in the simulation of the trigger and variation of the cuts used to select 
the data within the ranges allowed by the comparison between data and Monte Carlo 
simulations resulted in negligible changes in the corrected jet shapes. 

The systematic uncertainties for different values of r are correlated. For all the results the 
statistical errors are negligible compared to the systematic uncertainties. The total positive 
(negative) systematic uncertainty on ip(r) at each value of r was determined by adding in 
quadrature the positive (negative) deviations from the central value. The systematic uncer- 
tainties were then added in quadrature to the statistical errors and are shown as error bars in 
the figures. 



6 Inclusive jet photoproduction results 

6.1 Jet pseudorapidity dependence of the jet shape 

The jet shape for jets with Ej? > 14 GeV is presented in Figure [I] for four different rf et regions. 
It is observed that the jet shape broadens as if et increases. In the following the predictions of 
the PYTHIA generator are compared to the measurements. The predictions using HERWIG 
are very similar to those of PYTHIA and lead to the same conclusions. 



Comparison to parton shower model predictions 

The predictions of PYTHIA for direct, resolved, and direct plus resolved processes are com- 
pared to the measured jet shapes in Figure [I]. The admixture of the two processes was chosen 
according to the cross sections as given by PYTHIA. The dependence of the predictions on the 
specific sets of proton or photon parton distributions is negligible. The measured jet shapes in 
the range — 1 < if et < 1 are reasonably well described by the predictions of PYTHIA for either 
resolved plus direct or resolved processes alone. The predicted jet shapes in direct processes 
are significantly narrower than those of the data. These results are in agreement with the dom- 
inance of resolved processes in the E J T et range studied, as observed in the measurement of the 
inclusive jet differential cross sections [0, SI]. As rj^ et increases beyond rj^ et = 1, the predicted 



jet shapes increasingly deviate from those in the data, which are significantly broader in the 
most forward region (1.5 < rf et < 2). 

In order to study the relative importance of parton radiation and fragmentation effects, the 
predictions of PYTHIA for resolved plus direct processes without initial and final state parton 
radiation are also shown in Figure [I]. Since the predicted jet shapes based on fragmentation 
only are significantly narrower than the measured ones, it is concluded that, at the transverse 
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energies studied here, the shape of jets is mainly dictated by parton radiation and cannot be 
explained by hadronisation alone. 

Comparison to the parton shower model predictions of quark and gluon jets 

The difference between quark and gluon jets is modelled in PYTHIA by a leading-logarithm 
parton-shower approximation. The jet shapes, as predicted by PYTHIA, for quark and gluon 
jets are shown separately in Figure || together with the same data as in Figure |l|. The pre- 
dictions of PYTHIA for the two types of jets are computed separately for resolved and direct 
processes, and then averaged according to the cross sections given by PYTHIA. The predicted 
jet shapes are broader for gluon jets than for quark jets in each region of rf et and exhibit only 
a small dependence on nf et . 

The differences in the predictions for resolved and direct processes (see Figure [I]) come 
mainly from the different fractions of quark and gluon jets in the final state depending on 
the rf et region. This difference, in turn, originates from the different dominant two-body 
subprocesses. In the case of direct processes, where the dominant subprocess is photon-gluon 
fusion 7# — > qq, PYTHIA predicts the fraction of quark jets to be 80-100% depending on rf et . 
In the case of resolved processes, the dominant subprocess in the kinematic regime studied 
is q^g p — > qg and the predicted fraction of quark jets has a stronger T^-dependence: from 
~ 80% at 77 je * = — 1 to ~ 40% at rf et = 2. The comparison of measurement and prediction (see 
Figure shows that the broadening of the jet shapes in the data as rj^ et increases is consistent 
with an increasing fraction of gluon jets. Therefore, one possible reason for the deviation 
between data and prediction is an underestimation of the fraction of gluon jets in the region 
rf et > 1 by PYTHIA. 



Comparison to model predictions including multiparton interactions 

An excess of transverse energy outside of the jet cone for jets with rf et > 1 with respect to 



the expectations of PYTHIA was observed in previous studies || 31]. Since this excess may 



have some effect on the comparisons discussed so far, its influence on the jet shape has been 
considered. In order to simulate an increased energy flow the PYTHIA generator including 
MI was used; the cutoff (1.0 GeV) for the effective minimum transverse momentum for MI 
was tuned to reproduce the transverse energy flow outside of the jet cone in the data. The jet 
shapes predicted by PYTHIA MI are also shown in Figure |2|. Comparing with Figure [I], it is 
observed that the effects of MI on the jet shape are very small in the region —1 < r]^ et < 1, 
increase gradually with rf et and yield an improved description of the data in the region if et > 1. 
The differences in the predicted jet shapes between PYTHIA and PYTHIA MI are considered 
as an estimate of the uncertainty in the simulation of the energy flow outside the jet cone due 
to a possible underlying event. The conclusions that the measured jet shapes are consistent 
with the dominance of resolved processes and that the parton radiation is the main mechanism 
responsible for the jet shape still hold when this uncertainty is taken into account. 



Comparison to model predictions for fixed r 

Figure [3] shows the measured jet shape at a fixed value of r = 0.5, ip{r = 0.5), as a function 
of rf et . The predictions of PYTHIA for quark jets with and without MI span a band above the 
data, while those for gluon jets with and without MI span a band typically below the data. It is 
also seen that the effect of MI is larger as if et increases. The predictions for resolved plus direct 
processes are also shown in Figure |3|. The prediction of PYTHIA fails to describe the relatively 
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strong broadening of the measured jet shape for ry je * > 1. As mentioned earlier, one reason 
might be that the fraction of gluon jets in the region if ei > 1 is underestimated. However, 
when the effects of a possible underlying event are taken into account, using PYTHIA MI, the 
measured broadening in the forward region is accounted for with the default fraction of gluon 
jets. Note that, in any case, as Tf et increases the measured jet shape changes from a value close 
to the upper band (quark jets) to a value within the lower band (gluon jets). It is concluded 
that the broadening of the measured jet shape as rf et increases is consistent with an increase 
of the fraction of gluon jets independent of the effects of a possible underlying event. 



6.2 dependence of the jet shape 

The E^ 1 dependence of the jet shape is presented in Figure |j. It is observed that the jets 
become narrower as Ej? increases. For Ej? > 17 GeV the predictions of PYTHIA for resolved 
or resolved plus direct processes reproduce the data reasonably well. In the lowest E^ region, 
differences between data and the predictions are observed. Again the predicted jet shapes 
for direct processes are narrower than for the data. PYTHIA including resolved plus direct 
processes, but without initial and final state parton radiation, predicts jet shapes which are 
too narrow in each region of Ej? . These comparisons show again that parton radiation is the 
dominant mechanism responsible for the jet shape in the whole range of Ej 1 studied. 

The measured jet shape at a fixed value of r = 0.5, ip(r = 0.5), shows a moderate increase 
with E^, as seen in Figure [|. Note that the jet shapes have been measured in ranges of E J T et 
and the data points in Figure || are located at the weighted mean in each Ej? range. The 
predictions for the dependence of the jet shape on E^ 1 in resolved processes reproduce the data 
except for the lowest E^ 1 data point. The predicted jet shape for direct processes is narrower 
than the data. The effects of a possible underlying event are estimated using the predictions of 
PYTHIA MI. The inclusion of these effects improves the description of the data in the lowest 
Eif 1 data point, but otherwise does not alter significantly the .E^-dependence of the jet shape. 



6.3 Comparison to fixed-order QCD calculations 



Lowest-non-trivial order QCD calculations of the jet shapes |]19[ [2(| are compared to our mea- 
surements in Figures |6] and [7|. These predictions include resolved and direct processes, and use 
the CTEQ4 |32|] (GRV-HO) proton (photon) parton densities. Since the jet shape is computed 
only to the lowest-non-trivial orderQ, 0(a s ), the predictions are subject to relatively large un- 
certainties due to the strong dependence on the renormalisation and factorisation scales. In 
the calculations shown here, these scales have been chosen equal to E^. Since the calculations 
include only up to three partons in the final state, not more than two partons can build up 
a jet. As a result, the overlapping and merging issues of the experimental jet algorithm are 



not reproduced in the theoretical calculation Wj\ B3I. An attempt was made to simulate these 



effects by introducing an ad- hoc Rsep parameter [|TTJ : two partons are not merged into a single 
jet if their separation in the r\ — (f) plane is more than Rsep- The calculations of the jet shape 
shown in Figures |6] and [7] have been made for various Rsep values^: for two fixed values of 

2 The lowest-non-trivial order contribution to the jet shape corresponds to the NLO matrix elements of the 
hard interaction. The NLO contribution to the jet shape is not available due to the lack of the relevant next- 
to-next-to-leading order matrix elements, and in any case, as pointed out in [j33| , it cannot be safely calculated 
for the iterative cone algorithm used here. 

3 Although the value Rsep = 1 was suggested |}5) for the comparison between the measurements of cross 
sections and theoretical calculations, other values can be used in order to match the measured jet shapes. 
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Rsep = 1-4 and 2.0, and for the value of Rsep which best reproduces the data |2(*| . 

The fixed-order QCD calculations with a common value of Rsep — 1-4 reproduce reasonably 
well the measured jet shapes in the region —1 < rf et < 1 and in the region E^ 1 > 17 GeV. In 
the forward region, rf et > 1, and in the lowest E 3 ^ 1 region, 14 < E^ < 17 GeV, the calculations 
with Rsep = 1-4 show significant deviations with respect to the data. The discrepancy is very 
similar to that observed between the predictions of PYTHIA without MI and the data. However, 
a satisfactory description of the data can be obtained by leaving Rsep as a free parameter (for 
each interval in r] 3et ). A crosscheck with the data which overlaid jets from different events 
showed that the minimum distance at which two jets are resolved as two distinct jets depends 
very little on rf Gt . Therefore, the variation of Rsep with 77 je * cannot be justified on these 
grounds. The use of different values of Rsep to describe the data in the forward and lowest 
E^ 1 regions may be mimicking the effects of QCD higher orders and of a possible underlying 
event, which at present are not included in the calculations. In addition, the comparison 
between the measured jet shapes, which are corrected to the hadron level, and the fixed-order 
QCD calculations at the parton level is subject to the uncertainty of hadronisation effects. 



7 Dijet photoproduction results 

The jet shapes have been measured for each of the two highest E 3 ^ 1 jets in the reaction ep — > 
jet + jet + X. Both jets are required to have E^ > 14 GeV and —1 < rf ei < 2. The 
measurements have been integrated over four non-overlapping regions in r] 3et of each jet. The 
results for ip(r) are presented in Figure |]and are compared to the predictions of PYTHIA. The 
jet shape broadens as rf et increases and is narrower than that of the inclusive jet sample (see 
Figure |1]). The predictions for resolved or resolved plus direct processes describe reasonably well 
the measured jet shapes in all regions of if et . The predicted jet shapes for direct processes are 
narrower than for the data in the region t] jet > 0, as expected from the dominance of resolved 
processes in that region. The comparison of the predicted jet shapes with the data in the region 
rjjet > j ^ oeg no ^- t ne discrepancies observed in the inclusive jet photoproduction study. 
This difference is attributed to a suppression of the effects of a possible underlying event in the 
case of dijet events since the requirement of two high-i^ 6 ' jets increases the contribution from 
direct processes and, for resolved processes, decreases the leftover energy for the remnants. The 
broadening of the measured jet shapes in dijet photoproduction as r) 3et increases is adequately 
reproduced by the predictions of PYTHIA. 

In dijet photoproduction the contributions of resolved and direct processes can be separated 
j5], |7| by using the variable: 

Enjei _ ' 
jet s I 1 

2yE e 

where the sum runs over the two jets of highest E^ and yE e is the initial photon energy. This 
variable represents the fraction of the photon's momentum participating in the production of 
the two highest E^ 1 jets. The LO direct and resolved processes largely populate different 
regions of x® BS , with the direct processes concentrated at high values. 

The results for ip(r) are presented in Figure |9| for both x® B smaller and larger than 0.75. It is 
observed that the measured jet shapes for x^ BS < 0.75 are broader than those for x® BS > 0.75. 
For the region of x° BS > 0.75, the jet shapes as predicted by PYTHIA including resolved plus 
direct processes describe well the data. The predictions of PYTHIA for the region x® BS < 0.75 
reproduce the data reasonably well though the latter are slightly broader. The inclusion of 



x OBS _ ^jets-^T ° ^ 
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the effects of a possible underlying event as modelled with PYTHIA MI leads to an improved 
description of the data in the region of x° BS < 0.75 and has a negligible contribution in the 
region of x° BS > 0.75 (not shown). 

8 Summary and conclusions 

Measurements of jet shapes in inclusive jet and dijet photoproduction in e + p collisions at 
y/s = 300 GeV using data collected by ZEUS in 1994 have been presented. The jet shapes refer 
to jets at the hadron level with a cone radius of one unit in the i] — (f) plane and are given in the 
kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 277 GeV. Jets with E 3 T et > 14 GeV 
and — 1 < rf et < 2 have been considered. The dependence of the jet shapes on E^ 1 and if et has 
been studied: the jet shape broadens as r\> et increases and narrows as Ejp increases. In dijet 
photoproduction, the jet shapes have been measured separately for two regions of x° BS , the 
fraction of the photon's momentum participating in the production of the two jets of highest 
EjT 1 , x° BS < 0.75 and x° BS > 0.75. These subsamples are dominated by resolved and direct 
processes, respectively. The jet shapes in the region x° BS < 0.75 are systematically broader 
than those in the region x° BS > 0.75. 

Leading-logarithm parton-shower Monte Carlo calculations of resolved and direct processes 
have been compared to the measured jet shapes. The predictions based on resolved or resolved 
plus direct processes describe reasonably well the measured jet shapes in the region —1 < 
rf et < 1 for inclusive jet photoproduction and in the full region of rf et for dijet production. The 
predictions including only direct processes are narrower than those measured in the data. The 
removal of initial and final state parton radiation in the Monte Carlo calculations gives rise to 
jet shapes which are too narrow compared to those of the data. The results are in agreement 
with the dominance of resolved processes and indicate that parton radiation is the dominant 
mechanism responsible for the jet shape in the E^ range studied. 

Fixed-order QCD calculations cannot reproduce the measured jet shapes over the full kine- 
matic range with a single value of the Rsep parameter. 

The observed broadening of the jet shape as rf et increases is consistent with an increase of 
the fraction of gluon jets independent of the effects of a possible underlying event. 
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ZEUS 1994 




Figure 1: The measured jet shapes corrected to the hadron level, ip(r), for jets in the range 
above 14 GeV in different if et regions. The error bars show the statistical and systematic errors 
added in quadrature. For comparison, the predictions of PYTHIA for resolved (dashed), direct 
(dot-dashed line), and resolved plus direct processes (solid line) are shown. The predictions 
of PYTHIA for resolved plus direct processes without initial and final state parton radiation 
(dotted line) are also included (labelled by 'Only Fragment.'). 
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Figure 2: The measured jet shapes corrected to the hadron level, ip(r), for jets in the 
range above 14 GeV in different 7]^ et regions. The error bars show the statistical and systematic 
errors added in quadrature. For comparison, the predictions of PYTH1A including resolved plus 
direct processes for quark (dot-dashed line) and gluon jets (solid line), and those of PYTHIA 
MI for resolved plus direct processes (dotted line) are shown. 
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Figure 3: The measured jet shape corrected to the hadron level at a fixed value of r = 0.5, 
ip(r = 0.5), as a function of rf et for jets with Ej? larger than 14 GeV. The error bars show the 
statistical and systematic errors added in quadrature. For comparison, various predictions of 
PYTHIA including resolved plus direct processes are shown: quark jets (thin dot-dashed line), 
gluon jets (thin dashed line) and all jets (thin solid line). The corresponding predictions of 
PYTHIA MI are displayed with thick lines. 



16 



ZEUS 1994 




Figure 4: The measured jet shapes corrected to the hadron level, 4>(r), for jets in the r) jet 
range between —1 and 2 in different E^ 1 regions. The error bars show the statistical and 
systematic errors added in quadrature. For comparison, the predictions of PYTHIA for resolved 
(dashed), direct (dot-dashed line), and resolved plus direct processes (solid line) are shown. The 
predictions of PYTHIA for resolved plus direct processes without initial and final state parton 
radiation (dotted line) are also included (labelled by 'Only Fragment.'). 
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Figure 5: The measured jet shape corrected to the hadron level at a fixed value of r = 0.5, 
ip(r = 0.5), as a function of E^ 1 for jets in the if et range between —1 and 2. The error bars show 
the statistical and systematic errors added in quadrature. For comparison, various predictions 
of PYTHIA are shown: resolved (thin dashed line), direct (thin dot-dashed line) and resolved 
plus direct processes (thin solid line). The predictions of PYTHIA MI for resolved and resolved 
plus direct processes are displayed with thick lines. 
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Figure 6: The measured jet shapes corrected to the hadron level, ip(r), for jets in the E^ 1 range 
above 14 GeV in different r]^ et regions. The error bars show the statistical and systematic errors 
added in quadrature. For comparison, the predictions for the jet shapes (solid, dashed and dot- 
dashed lines) based upon the fixed-order QCD calculations by M. Klasen and G. Kramer with 
various values of Rsep (see text) are shown. 
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Figure 7: The measured jet shapes corrected to the hadron level, ip(r), for jets in the rf et range 
between — 1 and 2 in different E^ 1 regions. The error bars show the statistical and systematic 
errors added in quadrature. For comparison, the predictions for the jet shapes (solid, dashed and 
dot-dashed lines) based upon the fixed-order QCD calculations by M. Klasen and G. Kramer 
with various values of Rsep (see text) are shown. 
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Figure 8: The measured jet shapes corrected to the hadron level, ip(r), for each of the two 
highest jets in dijet photoproduction. Both jets are required to have > 14 GeV. The 
measurements are shown for four different regions in if et . The error bars represent the statistical 
and systematic errors added in quadrature. For comparison, the predictions of PYTHIA for 
resolved (dashed), direct (dot-dashed line), and resolved plus direct processes (solid line) are 
shown. 
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Figure 9: The measured jet shapes corrected to the hadron level, ip(r), for each of the two 



highest E^ 1 jets in dijet photoproduction separated according to x^ Bb . Both jets are required 
to have E^ 1 > 14 GeV and —1 < rj^ et < 2. The error bars show the statistical and systematic 
errors added in quadrature. For comparison, various predictions of PYTHIA including resolved 
plus direct processes are shown: for the region x° BS < 0.75 (PYTHIA without MI as the solid 
line and PYTHIA with MI as the dot-dashed line) and for x° BS > 0.75 (PYTHIA without MI 
as the dashed line). 
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